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Abstract. We investigate the implications for the Planck 
Surveyor of the recent sub-millimetre number counts ob¬ 
tained using the SCUBA camera. Since it observes at the 
same frequency as one of the higher frequency science 
channels on Planck, SCUBA can provide constraints on 
the point-source contribution to the CMB angular power 
spectrum, which require no extrapolation in frequency. 
We have calculated the two-point function of these sub¬ 
millimetre sources, using a Poisson model normalized to 
the observed counts. While the current data are uncertain, 
under reasonable assumptions the point-source contribu¬ 
tion to the anisotropy is comparable to the noise in the 
353 GHz channel. The clustering of these sources is cur¬ 
rently unknown, however if they cluster like the z ~ 3 
Lyman-break galaxies their signal would be larger than 
the primary anisotropy signal on scales smaller than about 
10 arcminutes. We expect the intensity of these sources to 
decrease for wavelengths longward of 850/j,m. At the next 
lowest Planck frequency, 217 GHz, the contribution from 
both the clustered and Poisson terms are dramatically re¬ 
duced. Hence we do not expect these sources to seriously 
affect Planck’s main science goal, the determination of the 
primordial anisotropy power spectrum. Rather, the poten¬ 
tial determination of the distribution of sub-mm sources 
is a further piece of cosmology that Planck may be able 
to tackle. 
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1. Introduction 

The study of the anisotropy in the Cosmic Microwave 
Background (CMB) has the potential to teach us a great 
deal about the background cosmology in which we live, 
about the formation of structure and about the early uni¬ 
verse. Because of this promise, ESA selected the Planck 
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Surveyoi0 as the third Medium sized mission of its Hori¬ 
zon 2000 Scientific Program. With its wide range of fre¬ 
quencies, superb angular resolution, and high sensitivity, 
Planck has been hailed as the definitive CMB anisotropy 
experiment. 

While much of the cosmological information is ex¬ 
pected to come from the angular power spectrum of pri¬ 
mary anisotropies, the sky contains more than this simple 
imprint of the inhomogeneities at last scattering. In par¬ 
ticular the brightness fluctuations should be dominated on 
the smallest scales by point sources. One of the collateral 
science goals of Planck is to produce an all sky catalogue 
of such sources over a wide range of frequencies. Recent 
observations suggest that there may be many more bright 
sub-millimetre sources than previously expected, and it is 
the purpose of this paper to explore the impact of these 
findings on the Planck mission. While several authors (see 
below) have looked at the one-point statistics of these 
sub-mm galaxies, little has been done on the two-point 
function (or power spectrum) of these sources. It is this 
statistic which is most familiar to the CMB community, 
and which we concentrate on here. 


The most exciting recent observations in the sub- 
mm waveband have come from the new Submillime¬ 
ter Common-User Bolometer Array (SCUBA; Holland et 
al. 1999) on the James Clerk Maxwell Telescope. A combi¬ 
nation of the properties of the sky and the galaxies them¬ 
selves make the SCUBA 850 /rm filter the optimal one for 
cosmology. This waveband corresponds closely with the 
353 GHz channel of the High Frequency Instrument (HFI) 
of Planck. The central frequencies are almost identical, al¬ 
though the Planck bandwidth will be considerably larger 
{Ah’lv = 0.25 rather than the ~ 0.1 of the SCUBA 850 /im 
filter). SCUBA has now been used to make several deep in¬ 
tegrations which have detected distant sources at 850 /xm 
(Barger et al. 199U Bales et al. 1998, Holland et al. 1998 


Hughes et al. 1998, Small et al. 1997), and this has radi¬ 
cally altered our expectations for the importance of dusty 
galaxies at high redshift. 


^ http://astro.estec.esa.nl/Planck/ 
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A summary of the source count observations is pro¬ 
vided in Table |^. These new data thus provide us with di¬ 
rect measurements of the number density of bright sources 
(albeit currently only over small patches of the sky, and 
a limited range of flux) at a frequency directly relevant 
to the Planck science channels. We shall work primarily 
at 353 GHz, though we shall also extrapolate these counts 
to nearby frequencies using models of the spectral energy 
distribution of the sources. 

There have been several estimates of how many sub- 
mm sources Planck might be able to detect (e.g. Bersanelli 
et al. 19961) , as well as estimates of other one-point statis¬ 
tics, often referred to as ‘confusion noise’ (Condon 1974, 
Blain et al. 1998). As far as the fluctuations are con¬ 
cerned, several studies have already been carried out on 
the implications of point sources for the measurement of 
CMB anisotropies. Much of this work dealt more specifi¬ 
cally with radio galaxy contributions at low frequency (e.g. 
Tegmark & Efstathiou 1996), or dealt with far-IR sources. 


1991 


1997 


Scut 

N 

95% CL 

Area 

Density 

(mJy) 



(Sq. Deg.) 

10V(Sq. Deg.) 

2.0 

5 

(1.8,10.9) 

1.6 X 10"® 

O.Z_i_2 

3.0 

2 

(0.3, 6.4) 

2.5 X 10"® 

o.8tS;I 

2.8 

12 

(6.4,20.2) 

6.8 X 10"® 


4.0 

6 

(2.3,12.3) 

2.4 X 10"® 

2-5lo.9 

5.0 

5 

(1.8,10.9) 

6.0 X 10"® 

O.SlS.a 

8.0 

4 

(1.2, 9.4) 

4.0 X 10"^^ 

boigi 

Table 1. 

, SCUBA point 

source observations at 850 


but concentrated on the uniform background contribution 
and the corr elatio n functi on (e. g. Franceschini et al. 

Bond et al. 1991 , Wang 1991 , Gawiser & Smoot 
Blain et al. 1998). Very similar studies have also been 


(353 GHz). We list the flux level to which each field was 
searched (generally 3 times the rms level), the number 
of sources found, the Bayesian 95% confidence level on 
the mean counts coming from Poisson statistics, the area 
covered, the number density of sources (with ±lcr errors 
also from Poisson statistics). These numbers were taken 
from Hughes et al. (1998), Barger et al. (1998), Bales et 
al. (1998), Small et al. (1997), Chapman et al. (in prepa¬ 
ration), and Holland et al. (1998), respectively. 


carried o ut for the X-ray background (see Yamamoto & 
Sugiyama |1998 , and r eferen ces therein), the optical back¬ 
ground (see Vogeley 1998 , and ref erenc es therein), and 
more exotic backgrounds (e.g. Scott 1993 ). The mos t rele - 
vant work are the recen t pap ers by Toffolatti et al. ( 1998 ) 
and Guiderdoni et al. (1996) which contain many useful 
results, and estimates of confusion noise, although they 
were written before the new series of SCUBA measure¬ 
ments. Our approach here also differs from theirs in that 
our predictions are based on straightforward extrapola¬ 
tions from observable properties. We deliberately avoid 
any modelling of the complex galaxy evolution process. 

To calculate the impact of these new data on Planck it 
will be necessary to model the number density of sources 
as a function of flux, N{> S). We shall take care to ensure 
that our model does not overproduce, when integrated, the 
far -infra red backgroun d (FI B) light detected by (Puget et 
al. 1996, Fixsen et al. 1998). At 353GHz the background 
is approximately 0.14 ± 0.04MJysr“^, although this is 
not a directly measured quantity, so the real error bar at 
353 GHz may be larger. We construct our fiducial model 
so that the integrated light from the sources contribute 
essentially all of this background. 


2. Counts and the Far Infrared Background 

Gu ided by models of galax y form ation (e.g. T offolat ti et 
al. 1998| , Guiderdoni et al. 1998 , Blain et al. 1998 ), we 
model N{> Scut) as a double-power-law. 


N{> Scut) = No 


1-f 


( 1 ) 


(see also Borys et al. 1998), where for convenience we take 
So = 10 mjy. Such a parameterization will certainly not be 
valid for all values of S, but it will be adequate for our 
purposes. Matching to the general behaviour of successful 
models, and normalizing specifically to the Hubble Deep 
Field (HDF) counts, we obtain ‘fiducial’ values for these 
parameters of Nq = 4.05 x 10®sr“^, a = 0.8 and j3 = 1.8. 
This fit, plus the data of Table |^, are shown in Fig. 
Some models show the counts to steepen even more at the 
bright end, but this has little effect on any results, since 
it is already steep enough that any upper flux cut is not 
very important. We also show in Fig. an estimate for 
the number of pixels over the whole sky at 353GHz, to 
indicate where we expect one source per pixel. We have 
assumed an oversampling by a factor of 10 pixels per beam 
as an illustrative number. With our adopted source counts 
model, Planck will then have about one 20 mJy source per 
beam (which then sets the basic level of ‘confusion noise’). 

In practice Planck will only be able to detect individual 
sources with at best a flux cut of S' > 100 mJy using data 
from the 353 GHz channel alone. With extra information 
from the higher frequency channels (as well as information 
from other instruments, at least in some regions of the 
sky), it should be possible to remove all sources down to 
a few X 10 mJy. The SGUBA counts constrain the model 
at somewhat lower flux levels than these, however in our 
model the counts at the SCUBA flux levels contribute 
significantly to the IR background and hence the CMB 
fluctuations if the sources are clustered (see below). In any 
case the precise flux cut for Planck is not currently easy to 
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Soul (mJy) (mJy) 


Fig. 1. The number of objects brighter than flux Scut as 
a function of S'cut- The curve is equation (|^) with our 
fiducial parameters. The data are taken from Table 
with the second and third points offset for clarity. The 
error bars are ±lcr (68% Bayesian confidence region) er¬ 
rors based on Poisson statistics. The horizontal dashed 
line is an estimate of the number of pixels on the whole 
sky for the Planck 353 GHz channel, assuming 10 pixels 
per 4.5' beam. 


estimate, and so we have erred on the side of conservatism; 
if the flux cut ends up being higher than we are assuming 
here, then the fluctuations will only be larger. 

By describing the counts in this phenomenological way 
we avoid any direct modelling of galaxy formation, evolu¬ 
tion and spectral synthesis. Currently there are too many 
free parameters in these ‘semi-analytic’ models to yield a 
great deal of insight. Instead we prefer to use simple model 
fits to observables on the sky, which are motivated by the 
current data. Because of this we are considering only the 
two-dimensional distribution of objects on the sky, with 
no requirement on the radial distribution. 

The contribution of these sources to the FIB is just 
the total flux per unit solid angle, or 


rFIB 



abu 


( 2 ) 


which can be integrated by parts to yield 

■S'cut 

NdS, - N{S,^t) Scut (3) 

(a little care has to be taken with minus signs, since con¬ 
ventionally N = N{S>Scut))- The faint end limit for con¬ 
stant slope is just a{l — a)~^N{> Scut)Scut- We show 
in Fig. ^ the contribution to the integrated background 
light as a function of S'cut- Notice that the sources at the 
flux levels probed by SCUBA contribute significantly to 
the background. As we will see below, it is those clustered 


Fig. 2. The integrated flux as a function of the upper flux 
cut, with So = 10 mJy, a = 0.8, /3 = 1.8, and normalized to 
the HDF counts. The h atched region is the FIB detection 
of Fixsen et al. ( 1998 ). Note that in our fiducial model 
much of the integrated flux comes from sources near the 
flux levels detected by SCUBA. 



Faint End Slope (a) 


Fig. 3. The total integrated fltrx; as a function of the faint 
end slope, a of equation (|^, with the other parameters 
fixed at their fiducial values. The hatched region is the 
FIB detection of Fixsen et al. ([19981 ). 


sources which contribute most to the background that may 
be of greatest interest to us here. In Fig. ^ we show the 
contribution to the FIB, integrating to Scut = oO; i-e. the 
total background, as a function of the faint-end slope a. 
The FIB was first detected by Puget et al. ( 1996| ), and 
has recently been measured by Fixsen et al. ( 19981 ). Their 
value is shown in Figs. ||, ^ as the hatched region. In our 
fiducial model the sub-mm sources account for all of the 
FIB. 
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Let us acknowledge that the real situation may be more 
complicated. The counts may come from a number of sep¬ 
arate populations, and so of course there could be features 
in the actual curve. In addition there is the possibility that 
some more diffuse emission contributes to the FIB, and is 
not accounted for in these counts. Indeed there are some 
early indications (Hughes et al. 1998, Borys et al. 1998) 
that the counts may be flatter at the faint end than the 
form we have adopted. Again we have been conservative 
here; lower faint-end slopes would require higher overall 
normalization in order to match the background, imply¬ 
ing stronger fluctuations. 


3. Power spectrum for sources 


We shall quote our results in ‘temperature’ units as is 
usual in CMB anisotropy studies. For fluctuations about a 
mean, the conversion factor from temperature to intensity 
(or flux) is 


dBij _ 2fc / /cTcmb \" 

^ ) 
99.27 Jy sr-i 




dT 


{e^ - 1)2 


fi K 


(4) 


(e- - 1)2’ 

where B^, is the Planck function, k is Boltzmann’s con¬ 
stant, X = hv/kBTcMB = 1^/56.84 GHz is the ‘dimension¬ 
less frequency’ and 1 Jy = 10“26 yVni“2Hz“^. 

If we assume that the number of sources of a given 
flux is independent of the number at a different flux, and 
if the angular two-point function of the point-sources is 
w{9), then the angular power spectrum, C^, contributed 
by these sources is 

rScMt H W n 

d5. 


= 


r™' + 


(5) 


assuming that all sources with S>Scnt are removed. Here, 
= f S dN/dS dS is the background contributed by 
sources below Scut as before. Following the conventional 
notation (essentially introduced by Peebles 1973| ), is 
the Legendre transform of the correlation function C(9) 
produced by the sources and we is the Legendre transform 
of w(0): 




'(2£ + l)CePe(cos&) 


(9) = —Y^(2£+l)w£Pe(cos9), 


( 6 ) 


(7) 


with Pe{cos6) the Legendre polynomial of order £. The 
first term in equation (|5|) is the usual Poisson shot- 
noise term (see Peebles |198C §46, or Tegmark & Efs- 
tathiou 19961 ), the second is due to clustering, assuming 
that the clustering is independent of flux. 

Integrating by parts the Poisson term can be rewritten 


pScut 

Ceiiy) = 2 NS,dS,- N{Scut) Sl 
^0 


( 8 ) 



Fig. 4. The Poisson component of the angular power spec¬ 
trum of the point sources, in dimensionless units, for 
a range of flux cuts (i.e. removing all sources brighter 
than Scut)- Higher flux cuts give larger fluctuation lev¬ 
els. To compare with the level of primary anisotropy ex¬ 
pected, the prediction for a standard CDM spectrum is 
also shown, normalized to COBE. The thick solid line 
is the expected contribution to the power spectrum from 
noise in the 353 GHz channel of the Planck HFI. 


At low flux it becomes q :(2 - a) ^N{>Scut)Slc,u if the 
counts have slope a at the faint end. 

The shot-noise component of the angular power spec¬ 
trum of our fiducial model is shown in Fig. along with 
the primary GMB signal and the projected noise in the 
Planck 353 GHz channel. The Poisson fluctuations are cal¬ 
culated using our fiducial model and assuming that we cut 
out all sources brighter than 5'cut = 200,100, 50,10 mJy. 
Note that much of the fluctuation power comes from 
sources fainter than 50 mJy. The uncertainty in the nor¬ 
malization of these curves is directly proportional to our 
uncertainty in the counts, Nq, which we have normalized 
to the data of Table |^. Glearly TVq is uncertain to at least a 
factor of 2. However, increasing the normalization by this 
amount would overproduce the FIB unless the faint end 
slope is also modified. 

As an aside we mention that though the noise and 
Poisson component of the sources appear to have similar 
power spectra, they are nonetheless quite different enti¬ 
ties. The sources are on the sky, and thus contribute to 
the flux in every observation of that pixel, whereas the 
noise varies from observation to observation and by as¬ 
sumption is uncorrelated with the signal in the pixel ob¬ 
served. Given many observations of a given direction on 
the sky (as expected for Planck), the noise properties can 
be separated from the sky signal, even if they have the 
same power spectra. The estimated instrumental noise can 
be subtracted from the measured power spectrum, with 
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Fig. 5. As in Fig. showing a flux cut of 100 mJy only, 
and including the component due to clustering. We have 
modelled the clustering using equation (p^. 


any point source contribution being evident as an extra 
Cl = constant component. We therefore expect that anal¬ 
ysis of the Planck data set will include fitting for an excess 
white noise component, which would be most likely due 
to unclustered point sources. 

We now turn to the other contribution to equation (i- 
Unfortunately there is essentially no information about 
the clustering of the SCUBA sources at present. Hence 
the most conservative assumption would be no clustering 
at all, but that is obviously unreasonable. Hence, although 
we tried to be as conservative as possible when discussing 
the shot-noise power spectrum, for clustering we will just 
make some simple guess. If we approximate w{0) oc 9~^ 
then wi oc Assuming either that the sources cluster 

like galaxies today or as Lyman-break galaxies (Giavalisco 
et al. 1998|) at z ~ 3, we would expect /3 ~ 0.8-0.9. As 


an illustrative example we can assume that the sources 
which make up much of the FIB cluster like these Lyman- 
break galaxies at z ~ 3. We suspect that this may in fact 
be close to reality, since it is equivalent to assuming that 
the population is a highly biased one, collapsing early. 
On the other hand the SCUBA sources will span a wider 
redshift range than galaxies selected by the UV-dropout 
technique, thereby washing out the angular correlations to 
some extent. But, certainly the SCUBA-type sources are 
likely to be more highly clustered than IRAS galaxies at 
low redshift. Our example can be probably considered an 
optimistic one in terms of clustered power. If the objects 
are less biased than the Lyman-break galaxies (LBGs) by a 
factor of b~^ then one reduces wi, and hence the clustered 
contribution to Ci, by b~^. 


With our assumption w{9) ~ and 


Wi 


LEG 


~ 10 


-6 


I 

100 


- 1.1 


(9) 


On large angles w[&) is expected to drop below the power- 
law behaviour assumed in equation . The scale of non¬ 
linearity approximately marks this transition. If we as¬ 
sume the power spectrum is a power-law with index n then 
wi cx Thus we expect that on larger angular scales wi 
will flatten and gradually turn over to wi oc £. To take 
this into account we cause wi to become constant^ for 
£ ^ 100, where we have chosen this multipole because it 
is a round number, and not because we think it has any 
physical significance. Explicitly we approximate using 


1 

wi 


,,LBG 

^100 


„LBG 


( 10 ) 


We show the amplitude of the clustering signal in 
Fig. I Note that the contribution due to source clustering 
dominates over the Poisson term on the range of angular 
scales relevant to CMB anisotropies. In fact if the clus¬ 
tering proves to be this strong then Planck may be able 
to measure the power spectrum of the IR sources over a 
range of angular scales. We imagine that such a compo¬ 
nent will be included as one of the foreground templates 
to be fitted for in a full Planck analysis - the template 
would include frequency dependance, as well as a power 
spectrum which is white noise with perhaps one or two 
additional parameters to describe the clustering part. Ex¬ 
actly how to model this component may change as we 
learn more from SCUBA and other instruments. Eurther 
understanding of the clustering of these sources is clearly 
an important direction for future research. 

One other issue is the variance in these power spectra 
estimates. Assuming that the point sources are a Poisson 
sample, it is straightforward to estimate the cosmic vari¬ 
ance associated with this component of the angular power 
spectrum^ The variance is a sum of two terms, one due 
to the finite number of modes sampled by any given Ci 
and the other from the Poisson nature of the process. So 
we have 


Cl = (11) 

and i6Cif = + I St^^dS.. (12) 

The first term of equation is the usual result for Gaus¬ 
sian fluctuations, the second term is the extra variance 

^ Following our argument above there is no reason why wi 
might not drop to low-^, rather than go constant, but for 
i ^ 100 the CMB signal dominates anyway, so our precise 
assumptions are not important. If we cause we to approach £ 
immediately then the clustered signal is down from Fig. ^ by 
an order of magnitude at £ = 10. Thus we caution the reader 
that this simple model may overestimate the clustered signal 
at low-^. A more complete treatment would require knowledge 
of the redshift distribution and the evolution of the clustering 
of the sources. 

^ If the full power spectrum is a sum of uncorrelated terms, 
the variance of the total is just the sum of the variances. 
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Fig. 6. As in Fig. |^, except at 217 GHz assuming the fluxes 
scale as We have kept the flux cut at 100 mJy to iso¬ 
late the effect of changing frequency, however in princi¬ 
ple the higher frequency channels could be used to isolate 
sources with 217 GHz flux substantially less than 100 mJy. 
This would lower the power spectra even further. 

associated with the Poisson sampling. If only a fraction 
/sky of the sky is observed then the first term is increased 
by while the Poisson term is unchanged. For now it 
seems safe to assume that the uncertainty in our mod¬ 
elling of the sources (i.e. the normalization and shape of 
dN/dS) is larger than the estimate of equation so 

we can neglect the latter. 

We expect the intensity of these sources to decrease 
for wavelengths longward of 850/im. Assuming the flux 
decreases as similar to the slope of the FIB, we have 
calculated the contribution at the next lowest Planck fre¬ 
quency: 217 GHz. As expected the contribution is con¬ 
siderably lower, as shown in Fig. If we use informa¬ 
tion on the spectra of individual galaxies which have been 
detected by SCUBA, (e.g. Ivison et al. 1998, Hughes et 
al. 1998), then the derived slope may be steeper still, lead¬ 
ing to lower 217 GHz contributions. Obviously at even 
lower frequencies, the signal will be correspondingly re¬ 
duced. At the next higher frequency channel, 545 GHz, 
the sources are obviously brighter. We And that the Pois¬ 
son contribution with a 100 mJy cut dominates the noise 
by more than 2 orders of magnitude over the range 2 ^ 
t ^ 10^, with the clustered contribution potentially larger 
still. 

4. Conclusions 

We have investigated the implications of the recent 
SCUBA number counts for the Planck Surveyor. Since 
it observes at the same frequency as one of the main sci¬ 
ence c hannels on Planck, SCUBA can provide constraints 
on the point-source contribution to the CMB angular 


power spectrum which require no extrapolation in fre¬ 
quency. While previous authors have investigated mainly 
one-point statistics of several SCUBA fields, it is the two- 
point function which has the most impact on the CMB sci¬ 
ence which will be done with Planck. We have calculated 
the two-point function of point-sources, using a Poisson 
model normalized to the observed counts. While the cur¬ 
rent data are uncertain, under reasonable assumptions the 
point-source contribution to the anisotropy is comparable 
to the instrumental noise in the 353 GHz channel. We have 
emphasized that if the instrumental noise power spectrum 
can be accurately estimated (as is expected to be the case 
for Planck, since each pixel is observed multiple times), 
then even a white noise contribution from point sources 
could still be detected. The clustering of these sources is 
extremely uncertain, however if they cluster like the z ~ 3 
Lyman-break galaxies their signal would be larger than 
the primary anisotropy signal on scales smaller than about 
10 arcminutes. We expect the intensity of these sources to 
decrease for wavelengths longward of 850/rm. At the next 
lowest frequency channel, 217 GHz, the contribution from 
both the clustered and Poisson terms is dramatically re¬ 
duced. 

The bottom-line is that the sub-mm sources revealed 
by SGUBA will not have a strong impact on the most im¬ 
portant goal of the Planck mission, that of precisely char¬ 
acterising the CMB anisotropy. For the entire Low Fre¬ 
quency Instrument, and the three lowest frequency chan¬ 
nels of the HFI there will be no significant contribution. 
And certainly the signals in the higher frequency channels 
can be used to remove point sources and recover most of 
the CMB information even at 353 GHz. Moreover, the pos¬ 
sibility of actually measuring the Poisson and clustering 
signals over most of the sky for these galaxies provides 
Planck with yet another way of tackling fundamental cos¬ 
mological issues. 
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